CONSTRUCTING ARITHMETIC SUBGROUPS OF UNIPOTENT 

GROUPS 

WILLEM A. DE GRAAF AND ANDREA PAVAN 



Abstract. Let G be a unipotent algebraic subgroup of some GL m (C) defined 
over Q. We describe an algorithm for finding a finite set of generators of the 
subgroup G(Z) = G(l GL m (Z). This is based on a new proof of the result 
(in more general form due to Borel and Harish-Chandra) that such a finite 
generating set exists. 

1. Introduction 

Let G be an algebraic subgroup of GL m (C) denned over Q, where m > 1. Then 
for a subring R of C we set 

G(R) = Gn GL m (R)- 

The group G(Z) and any other subgroup T of G(Q) commensurable with it are 
called arithmetic subgroups of G. 

Arithmetic groups occur in many contexts. Examples are: the automorphism 
group of a finitely generated nilpotent group ([10]. Chapter 6), the group of units 
of the ring of integers of a number field, and the group of units of the group algebra 
ZQ, where Q is a finite group. A celebrated theorem of Borel and Harish-Chandra 
([3]) says that arithmetic groups are finitely generated. In this paper we consider 
the problem of computing a finite set of generators of an arithmetic subgroup of a 
unipotent group. 

This problem was also treated in the paper [5] by Grunewald and Segal, where 
a general algorithm for all arithmetic groups was outlined. However, their declared 
aim was to show that such a computation is, at least in principle, feasible, and no at- 
tempt was made to make the algorithms as efficient as possible. And unfortunately 
it appears to be extremely hard to use their algorithm in practice. 

In this paper we describe a practical algorithm for finding a finite set of gener- 
ators of G(Z), in case G is unipotent, that is to say, all its elements are unipotent 
matrices. As a byproduct this yields an independent proof of the Borel-Harish- 
Chandra theorem in this case. Also, we can show that the groups Gl are T-groups 
of Hirsh length equal to dim G. In order to show that the algorithm is practical 
we have implemented it in the language of the computer algebra system GAP4 (cf. 

!)• 

We now sketch the main idea of the algorithm. Let V be the vector space on 
which G acts naturally. Let 

= V < Vi < ■ ■ ■ < V n = V 

be a flag of V with respect to the action of G (this means that for v £ Vi and g G G 
we have gv = v mod Vi-%). Then we can form the G- module V* — V n -\ © pj-- In 
informal terms the matrix of a g G G acting on V* is formed from the matrix of 
its action on V by taking the block in the upper left part of the matrix, and the 
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block in the bottom right part of the matrix, and constructing the block matrix 
consisting of these two blocks. Now let Q be the image of G in GL(V*); then we can 
recursively compute generators of Q(Z). The recursion works because the Q-flag in 
V* has smaller length. Let tt : G — > Q be the projection. In Section [2] we describe 
7r(G(Z)) (Proposition ^. 7|l . We show how to find generators of 7r(G(Z)), and their 
preimages in G(Z). Let N(Z) C G(Z) denote the kernel of ir. We find a finite set 
of generators of -/V(Z), and joined to the elements of G(Z) found earlier this solves 
the problem. 

These ideas are detailed in Section [2l In Section [3] we illustrate them with a 
simple example. The constructions of Section [5] do not immediately yield an implc- 
mentable algorithm. In order to obtain that we need some technical preparation. 
In Section U] we describe some results that allow us to work with the Lie algebra 
of G rather than with G itself. Section [5] contains some material on T- groups. In 
Section [6] we describe some algorithms for lattices that we need. Then in Section 
[7] we give a detailed description of the main algorithm, and prove its correctness. 
The last section describes some practical experiences with our implementation of 
this algorithm in GAP4. 

2. The derived representation 

The goal of this section is to introduce some notation, and to prove the results 
that underpin the main algorithm. 

Let V be a finite dimensional vector space over Q, and L a full-dimensional 
lattice of V. We do not prove the next lemma here; it will follow from Lemma [ 



Lemma 2.1. Let U and W be two subspaces ofV with U C W. Then there exist 
subspaces U' and W of V such that 

W C U' 

and equalities 

U®U' = V = W@W 

and 

(l nu) + (l rw') = l = (l nw) + (l n W) 

hold. 

In the notations of the lemma above, we say that W C U' is a system of L- 
complements for U C W . 
Now let 

= V < Vi < ■ ■ ■ < V n = V 
be a chain of subspaces of V with n > 1. Then we consider the vector space 

V* = V n -i © 

We call it the derived vector space. Also, we have the full-dimensional lattice 

vi 

of V* , which we call the derived lattice, and the chain of subspaces 

= v$ < vs* < • • • < v*_i = V* 

of V* where 

V* = Vi-i 



v l+1 



which we call the derived chain. Note that its length is n — 1, which is strictly less 
than the length of the chain of V. 
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Now let W n -i C W\ be a system of L-complements to Vi C V^_i, and let us 
denote by V V\ the projection of V onto Vi along W\. By Lin(W n _i, Vi) we 
denote the space of all linear maps W n —i ~ > Vi. Then we consider the map 

(1) e : End(V) — ► Lin(W n _i, Vi) ip h-> 7r o <p\w n _ x - 

We refer to it as the error map induced by the system W n -\ C Wi. Further we 
define 

(2) r = {7 e Lin(W„_i, Vi) 1 7 (L n Wn-i) c L n Vi}. 

Since L n W„_i is full-dimensional in W n -i, T is a lattice of Lin(W n _i, Vi), and 
it is full-dimensional since L n Vi is full-dimensional in Vi . We refer to it as the 
lattice induced by the system W n -i Q W\. 

Now let G be a unipotent algebraic group defined over Q acting faithfully on V, 
and suppose that 

= V Q < Vi < • • • < v n = v 

is a flag of V with respect to G, i.e., for all v £ Vi we have gv = v mod Vi_i for all 
g E G. We consider the subgroup 

G L - {9 6 G(Q) | 5 L = L} 

of G(Q). We want to find a finite set of generators of this group. 

Since Vi and V^_i are G-stable subspaces of V, G acts on both V n -i and 
hence on their direct sum, that is to say, on the derived vector space. We refer to 
the action of G on V* as the derived action. Further, we denote by N its kernel, 
which is of course a unipotent algebraic group over Q acting faithfully on V , by Q 
its image, which is a unipotent algebraic group over Q acting faithfully on V* , and 
by 7r the projection of G onto Q. The following lemma is well-known; it follows 
directly from the commutativity of diagram ((3|) in Section [4j 

Lemma 2.2. The projection n maps G(Q) surjectively onto Q{Q). 

Also we set 

N L = {g E N(Q) I gL = L) 

and 

Ql* ={qeQ(Q)\qL* = L*}. 
Of course the derived chain is a flag for V* with respect to the action of Q. Since G 
acts faithfully on V we can regard elements in G(Q) as automorphisms of V . The 
same consideration applies to elements of N(Q). So we can apply the map e to the 
elements of these groups. 

Proposition 2.3. For every g E G(Q) and h E N(Q) we have 

e(g ■ h) = e(g) + e(h), 

where e is as in (1$). 

Proof. Let v E W n -\. Since h is an automorphism of V acting as the identity on 
v_ 

Vi ' 

h(v) - v G Vi.. 

Further, g is an automorphism of V acting as the identity on Vi, hence 

g{h(v) — v) = h{v) — v. 

Since W n -\ C Wi, we have ir(v) = 0. Hence applying ir to both sides of the 
previous identity and using linearity we obtain 

7T o g o h{v) = tt o g(v) + ir o h(v) 

hence the thesis. □ 
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Of course, N(Q) acts on G(Q) by multiplication on the right. Once we endow 
Lin(W n _i,Vi) with the obvious group structure given by addition, the previous 
proposition implies that the restriction of e to N(Q) is a group morphism. Hence 
N(Q) acts on Lin(W„_i, V x ) by 

LmiWn-uVx) x N{Q) ^hiniWn-uVi) (x, h) » x + e(h). 

With these observations, we can restate the previous proposition saying that the 
restriction of e to G(Q) is a morphism of A r (Q)-sets. 

Now let us denote by W the image of N(Q) through e. Then Proposition 12.31 
gives us the 

Corollary 2.4. There exists a unique map e : Q{Q) — ► Lin(wV-i,Vi) swc ^ ^ g 
diagram 

— £Lra(W n _i,Vi) 



Lin(yy„-i,yi) 
w 



is commutative. 



Proof. By Lemma |2~2"| the map ir : G(Q) — > Q(Q) is surjective. Also, if g, g' € 
are such that Tr(g) = Tr(g') then g^ x g' G N(Q), hence due to Proposition 12.31 we 
obtain 

e(g') = c(ffs-V) = <9) + eQrV) 

thus 

e(g) +W = e(g') + W. 
These two facts show that the function 

e : Q(Q) — > Un (Wn^i,Vi) q ^ t{g) + w 

where g is any element of G(Q) such that ir(g) = q, is well defined and, of course, 
it makes the diagram above commutative. If e' is another such a function, then 

i O TT = i' O TT 

hence, by surjectivity of tt it follows that e' = e. □ 
Also we set 

G L * = {geG(Q)\gL* = L*}, 
In other words, an element g G G(Q) lies in Gl* if and only if ir(g) lies in Ql* . Of 
course, Gl* contains both Gl and N(Q). 

Lemma 2.5. Let g G Gl* ■ Then g G Gl if and only if e(g) G T. 
Proof. If g G G L then g(L n W n _i) C L, hence 

vr o 5 (l n W n -i) c ?r(L) = 7r((i n Vi) © (L n Wi)) = inVi 

hence e(<?) G T. 

Now let g G Gl* such that e(<?) G T. Then g is an automorphism of V fixing 
both L + V\ and L n V^_i, and such that nog sends i n W n _i in L n Vi. In 
particular, since the preimage of £ D Vi through 7r is W\ + (LdVi), we have that 
g(LnW n -i) C Wi + (£nVi). Further, since LfW n _i C L + Vi and 5 fixes L + 
we have that g(L n W„_i) C L + Vi. Hence 

<?(l n w n -i) c (l + Vi) n (Wi + (i n Pi)) 

Since L = (LnVi) + (Lfl W\), applying Dedekind's modular law we obtain equality 

L= (l + Vi) n (Wi + (in Vi)) 
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hence g(L n W n -i) C L. Since L = [L n K-i) + (in W„_i) and 5 fixes L n V„-i, 
this shows that g(L) C L. Now let Z £ L. Since 5 fixes V\ + L, there exist l\ G Vi 
and I2 € L such that <?(/i) + .9(^2) = Z. Since g(£) C L, we have that g(h) G £) 
hence in particular that g(Zi) = Z — fj(h) G £• Since 5 fixes Vi, we also have that 
G Vi, hence g(h) E L n Vi. Since Vi C K-i, also gr(/j) £ Lfl K_i holds. 
Since g fixes L Pi V n -i, we obtain that Zi G L n V^_i, hence Zi + Z2 G £, hence 
L C g(L). So 5 G G L . □ 

Proposition 2.6. The map given by the chain 

G L * — ► Lm(VF n _i, Vi) — > 

is a group morphism with kernel Gl ■ 



Proof. Let /, g E Gl* ■ Then they are automorphisms of V acting as the identity 
on Vi, fixing V n -i and acting as the identity on v v - . Further, they fix L n V n -i 
and L + Vi. Now let Z E L. Of course, g{l) - I E K-i Pi(L + Vi). Since Vi C V n -\, 
applying Dedekind's modular law we obtain V n -i fl (L + Vi) = Vi + (L n 
which shows that 

f(g(l)-l)-(g(l)-l) GLHVn-! 

hence 

f°g(i)-f(i)-g(i)£L. 

Since 7r(L) = L n Vi, we finally obtain 

7T o f o g(l) - 7T o /(Z) - 7T o g(Z) G i n Vi 

which shows that the map is a group morphism. By Lemma 12.51 its kernel is 
G L . □ 



Proposition 2.7. Let e be as in Corollary \2.4\ The map ^ given by the chain 

t LinfW^Vi) Lin(W n -i,yi) 
y W ' W + T 

is a group morphism. Its kernel is equal to the image of Gl through n. 

Proof. Since ir : G(Q) — » Q(Q) is surjective and g E G(Q) is in Gl* if and only if 
^(9) G Ql*, we obtain a surjective map 7r : Gl* — * Ql*- By commutativity of the 
diagram in Corollary 1 2. A\ we also obtain the commutative diagram 

G L * Lin(W„_i, Vi) iW-i.Vi) 



By Proposition I2.6[ the top row is a group morphism, hence also the bottom row 
is. Again by Proposition 12. 6[ Gl is the kernel of the top row. Hence, since the 
diagram above is commutative, the image of Gl through tt lies in the kernel of \&. 
Now let q be in the kernel of 'J, and let g E Gl* be a preimage of q through tt. By 
commutativity of the diagram above, we have 

e(g) ET + W 

Now let w G W such that e(g) + w E T. Since W is the image of N(Q) through e, 
there exists h E N(Q) such that e(h) — w. Hence by Proposition 12.31 we have 

e(g-h)=e(g) + e(h) GT 

thus by Lemma l2~5l we obtain that g ■ h E Gl- Of course, n(g ■ h) = q. □ 
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3. An example 



Let us consider 



G 



( 



\ 



2° 
C 

1 



G GL 4 (C) such that a,b,c e 



/ 



It is easy to check that G is an algebraic subgroup of GL 4 (C) defined over Q. Since 
it is contained in the set of upper-unitriangular matrices of GLi4(C), it is unipotcnt. 
In this example it is rather straightforward to find a set of generators of G(Z) 
directly. However, in this section we illustrate the results of the previous section 
by showing how they help us finding a finite set of generators for G(Z). 

G acts faithfully on Q 4 by matrix- vector multiplication; also, Z 4 is a full-dimensional 
lattice of Q 4 . Thus we can consider the subgroup G Z 4 of G, and it is easily seen 
that 

G(Z) =G Z 4. 

The chain of subspaces 

= V < (ei, e 2 ) = Vx < (e u e 2 , e 3 ) = V 2 <V 3 = Q 4 

where e\, e 2 , e 3 and e 4 are the standard basis of Q 4 , is a flag of V with respect to 
the action of G. The derived vector space has basis given by (ei, 0), (e 2l 0), (e 3 , 0), 
(0, e 3 + V\) and (0, e 4 + Vi); through this basis we can identify it with Q 5 , and 
under this identification the derived lattice corresponds to Z 5 . The kernel of the 
derived action is 

/ 1 b \ 

G GL 4 (C) such that b e C 



N 



In particular, 



e GL 4 (C) such that b e Z 



and it is straighforward to check that it is an infinite cyclic group with generator 



The subgroup 



Q 
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€ GL 5 (C) such that a, c € C > 
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c 
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of GL 5 (C) is the image of the derived action, the projection of G onto Q being 



/ 



:G^Q 



( 
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Now Qjfi is a torsion free abelian group of rank 2 with basis given by 
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1 J 




V o 











1 / 



A system of Z 4 -complements for V\ C V2 is given by W 2 C Wi where 

W 2 = ( ei ) Wx = {e z ,e 4 ) 

Using the basis for W 2 and the basis e\, e 2 for V\ we can identify Lin(W2, V\) with 
M2xi(Q); under this identification, the induced lattice T corresponds to M2 X i(Z). 



Further, using the standard basis of ( 
this way, the error map is 



we can identify End(Q 4 ) with M 



4xM 



In 



M 4x4 (Q) ^M 2x i(Q) 



«1,4 
02,4 



ai.i 

02,1 
03,1 
V fl 4T 

The image of the rational points of N through e is the subspace of M2 X i(Q) gener- 
ated by (J). So in the notation of Section[2]we have W = ((q)) and V = M2 X i(Z). 
Furthermore, the map 4" from Proposition 12.71 goes from the rational points of Q 



ai,2 
02,2 
03,2 
04,2 



«1,3 
d2,3 

03,3 

^4,3 



Ol,4 
02,4 
03,4 
04,4 



to 



M 2 



W+T 



Now we need two matrices g\ and g 2 in G(Q) whose images through 7r are q\ 
and q 2 , respectively. Their existence is guaranteed by the surjectivity of 7r. For 
example, we can take 



9i = 



1 



V 



92 



t 



\ 



This shows in particular that g\ and g 2 are in G-i 
the diagram in Corollarv l2.4( we have that 



Also, using commutativity of 



+ W + T 



*(«u) 



W 



Now the kernel of ^> is generated by qi and q\. Therefore, by Proposition 12.71 
tt{Gh) is generated by gi and q\. Their preimages, g\ and g\ are only guaranteed 
to lie in Gjp \ however, here we see that they are alreay in Gji. Now the kernel of 
7r restricted to G Z 4 is 7V Z 4. So giN Z 4 and g 2 N^4 generate jfj. We conclude that 

n,gi,gl generate G Z 4. 

Of course we could have made a different choice for the preimages of q\ and q 2 . 
For example, we could have taken 



f/i 



instead of g\ . Then g[ is in Gi 
kernel of we get e(qi) € W 
Corollary e(g[) S W 



1 



V 

but it is not in G Z 4. However, since gi is in the 
- r and by the commutativity of the diagram in 



r. This can of course also be checked directly as 
<9'i) - 
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Now we note that 



/ 1 i \ 

10 

10 

y o o o i ) 



with e(n') — ( ^ ^ 



The existence of such an n' is guaranteed by the fact that W is the image of the 
rational points of N through e. Since N is in the kernel of it we have Tr(g[ ■ (n')^ 1 ) = 
ir(g'i) — qi- So we can work with <?5.(n') _1 as preimage of q\. This choice of preimage 
works for us as, due to Proposition 12.31 we have that 



e(^.(n')- 1 )=e( 5l )-e(n')=( o)~(o) = (o) Gr ' 

Hence g'^n 1 )^ 1 G G%4 and we can apply the same considerations as above to prove 
that n, g'^n')^ 1 and g\ are a generating set for Gg4. As a matter of coincidence, 
we note that 

Finally we note that n commutes with both g\ and g\ , hence the chain of subgroups 

1 < (n) < (n,g-i) < (n,gi,#f) = G Z 4 
is a central series for G Z 4 with infinite cyclic factors. 

4. The Lie algebra connection 

Let G C GL m (C) be a unipotent algebraic group. As illustrated in Section [3l 
the results in Section [2J in principle yield an algorithm for finding a finite set of 
generators of G(Z). However, to make it work efficiently in practice we rather work 
with the Lie algebra of G than with G itself. In this section we describe the main 
reults that we need for that. 

First we review some standard facts on the Lie algebra of an algebraic group; 
for more details we refer to [2J, [5], [12 | .[13 | . 

As customary we denote the Lie algebras of the algebraic groups G,H, . . . by 
0, fj, . . .. Let G be a unipotent algebraic group defined over Q, acting on a vector 
space V. Then also acts on V. Furthermore, G is connected, and hence a subspace 
U C V is G-stable if and only if it is 0-stable. If this is the case then we get a 
G-action and a 0-action on U, and those are compatible, in the sense that the 
corresponding 0-representation is the differential of the G-representation. Similarly 
we get compatible G- and 0-actions on quotients and direct sums of modules. 

An important role in our algorithm is played by the exponential mapping. For 
a nilpotent x G g! m (C) we set 



exp(a;) 

and for a unipotent u G GL m (C) 



n— 1 ^ 

EX 

i=0 



1 



iog(t*) = 2^(-i) 



Since G is defined over Q we have that C 0l m (C) has a basis such that all 
elements have coefficients in Q. The Q-span of such a basis is denoted 0q. Then it 
is well-known that the maps exp : 0q — » G(Q) and log : G(Q) — ► 0q are mutually 
inverse. In particular, when we work with the Lie algebra 0q we keep control over 
the elements of G(Q) by these mappings. 
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Now let W be another finite dimensional vector space over Q, H a unipotent 
algebraic subgroup of GL(W), and ip : G — > a morphism of algebraic groups. 
Then we have the diagram 



1. 'U 



(3) 



cxp 
log 



dip 



H( 



cxp 



and it turns out that it is commutative, i.e., exp(d(p(x)) = (p(exp(x)) for all x G g 
(cf. [5], Chapter V, §4, Proposition 15). 

Now we return to the setting of Section [2] A sequence of subspaces 

= V < Vi < ■ ■ ■ < V n = V 

is a flag for the action of G if and only if it is a flag for the action of 0. (The 
latter means that g-Vi C Vi~\ for i > 0.) In particular, acts on the derived 
vector space V*, and the corresponding representation of g is the differential of the 
representation of G onV*. In particular this means that n, which is the Lie algebra 
of N, is the kernel of d7r and q = d7r(fl), where q is the Lie algebra of Q. 

Proposition 4.1. n is central in g. 

Proof. Let x G g and y G n. Then x is an endomorphism of V such that x.V G V n _i 
and x.Vi — 0, and y is an endomorphism of V such that y.V C V\ and y.V n -\ = 0. 
Now let v £ V. Then y.v G V\, hence x.y.v = 0. Also, y.v G Vi, hence x.y.v = 0. 
Thus [x, y].v = 0. Since g acts faithfully on V, [x, y] = 0. □ 

Now let L be a full-dimensional lattice of V, W n -i Q W\ a system of L- 
complements to V\ C V n -i, and let us denote by e and by T the induced error 
map and the induced lattice, respectively. Since g acts faithfully on V, we can 
regard its elements as endomorphisms of V. The same consideration applies to n. 
So we can consider the restriction of the map e to and n. 

Proposition 4.2. The restriction of the induced map e to n is injective. 

Proof. Let x G n such that e(x) = 0. Since x(V) C Vi, for every v G V we have 

= 7r o 



Since e(x) — 0, n o x(W n -i) 
finally have x(V) = 0. 



0. Hence x{W n -\) = 0. Since x{V n -i) = 0, we 

□ 



Further, we define 

n L = {x G n | e(x) G F}. 

By the previous proposition, it is a full-dimensional lattice of n. Since n acts 
faithfully on V and V admits a flag with respect to this action, then n, regarded 
as a Lie subalgebra of gl{V), consists of nilpotent endomorphisms. Hence we can 
consider the diagram 

° xp GL(V) 




Lm(W n _i,Vi) 



Proposition 4.3. The diagram above is commutative. 
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Proof. Let us denote by idy the identity endomorphism of V. Since W n -\ C W\, 
e(idy) = 0. Now let x G n. Then x(V) C F x and x{V n -x) = 0, hence x 2 (V) = 0, 
and 

exp x = idy + cc 

Thus 

e(exp x) = e(idy + ar) = e(idy) + e(x) = e(x) 

□ 

5. T-GROUPS 

The groups G(Z) that we are after are finitely-generated nilpotent and torsion 
free. Such groups are called T-groups in the literature (cf. [10]). In this section we 
review some facts that we need on T-groups. 

It is known that any T-group admits a (proper normal) central series with infinite 
cyclic factors; the other way round, every group admitting such a series is clearly 
a T-group. Now let G be a group, and let gi, . . . , g n be a (ordered) set of elements 
of G. Then we can consider the chain of subgroups 

G\ > G2 > • • • > G n > G„+i = 1 

of G where for every i = 1, . . . , n, 

Gi = (gi, ...,g n ) 

We call it the chain associated to g\,...,g n - Further, we say that gi, ■ ■ ■ ,g n is 
a T-sequence for G if the associated chain is a proper central series for G with 
infinite cyclic factors. It will be convenient to extend this terminology saying that 
the empty set is a T-sequence for the trivial group. Every T-group G has a T- 
sequence, and the length of a T-sequence is an invariant of the group, called the 
Hirsch-length of G. Furthermore, if gi, ...,g n is a T-sequence for G then every 
element g G G can be written g = g^ 1 ■ ■ ■ gf™ , where ej G Z. 

Now let A be an abelian group, and let a%, . . . , a n G A, Then we consider 

L = {(ei,...,e„) GZ"|a^...<" = 1} 

which is of course a subgroup of Z". We call it the relation lattice of ai, . . . , a n in 
A. For every e = (ei, . . . , e„) G L, e ^ 0, we can define its height as the minimum 
j = 1, ... ,n such that ej ^ 0, and its leading coefficient as the integer ej. Now let 
e^ x \ . . . , e 1 ™) be a basis of L, where for j — 1, . . . , m 

e W = ( e «,...,eW). 
We say that the basis is in Hermite normal form if the matrix 

/ ■ • • e« \ 



G M mx „(Z) 



is. This means that there exists a sequence of integers 

!<«!<•••< i m < n 



e ?' } = 



such that 

for all j = 1, . . . , to and all 1 < i < z 3 , and that 

< ej } < eg 5 

for every 1 < k < j < m. It is known that L admits a unique basis in Hermite 
normal form. We note that the height of any non-zero element of L is any of the 
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integers i\, . . . , i m . Further, if its height is ij for some j = 1, . . . , m, then its leading 
coefficient is a (non-zero) multiple of ef' . 

The next lemma is a somewhat stronger version of a result of Eick (cf. [6], 
Lemma 3.19). 

Lemma 5.1. Let G be a T -group, A an abelian group, and let tp : G —* A be a mor- 
phim of groups. Further, let g\, . . . , g n be a T -sequence for G, and let . . . , e^™ 1 ' 
be the basis in Hermite normal form of the relation lattice oftp(gi), . . . , tp{g n ) in A. 
For j = 1, . . . , m set 

kj = Si 1 ■■■9n ■ 
Then k±, . . . , k m is a T -sequence for the kernel of tp. 

Proof. Let us denote by 

K x > K 2 > ■ ■ ■ > K m > K m+1 = 1 

the chain of subgroups of G associated to fci , . . . , k m . We want to show that it is a 
(proper normal) central series for ker(/3 with infinite cyclic factors. Since the basis 
e^ x \ . . . , e( m ) is in Hermite normal form, we have the sequence 

1 < i\ < i2 < ■ ■ ■ < im < n 

defined as above. It is convenient to set io — and i TO +i = n + 1. Also, let us 
denote by 

G = G\ > G-2 > • • • > G n > G n +i = 1 

the proper central sequence with infinite cyclic factors for G associated to gi, . . . , g n . 
Then it is enough to prove that for every j = l,...,m + lwe have 



kertp n Gi 3 _ 1+ i = ■ ■ ■ = kevp (~l Gi 3 -i = ker^s n G i} = K 



j 

Indeed, suppose that the previous equalities hold. Then A'i = kerip n G^+i = 
keiip n Gi = kenp D G = kenp. Also, for every j = 1, . . . , m+ 1, Gi-<G, hence Kj — 
keiip Pi Gi . < kertp n G = kertp. This shows that the chain is a normal series for kertp. 
Further, for every j — 1, . . . , m, the map kertp — > — has kernel kertp n Gij+i — 

Kj+i- Hence it factors through a group monomorphism j? 1 ^ — > G G +i . The image 

t Kj u .. • + 1 (ker^nG* )G« +i (kor^G, +1 )nG,/ G, 

of g 3 through it is — ^ — 1 — = ^-^ 1 — = pr 2 < ^ , and 

"3 + 1 ' trij + l Lri i + 1 ~ Lri j + 1 

the image of kjKj + ± is g i i Gj.+i. This shows that the series is central and with 
infinite cyclic factors. 

So we have to prove the previous equalities. It is clear that for every j = 
1, . . . , m + 1, 

kertp n G i3 _ 1+ i D ■ ■ ■ D kertp n G^-i 2 kertp n G i] D Kj 

and it remains to prove the reverse inclusions. We proceed by induction on j. 
Let us consider the base case j = m + 1. Then Kj is trivial, and all we have to 
show is that for every I = i m + 1, . . . , n + 1, kertp n Gi is trivial, too. Again, we 
proceed by induction on I. In the base case I = n + 1, it is obviously true. Now let 
J = im + 1, . . . j 71 and suppose that kertp n Gj+i is trivial. Let <? G kertp n G;. Then 
g = gfh for some e € Z and ft, G G/+i, and 

tp( ff; ) e + ^W -o 

Since ft. G {gi+i, ■ ■ ■ , gn), then tp(ft) G (<p(<?;+i), ■ ■ ■ , <p{gn))- Thus if e 7^ 0, then there 
would exist an clement in the relation lattice with height I, which is impossible. 
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Hence e = 0, hence g = h G kenp n G;+i, hence g = 1 by the inductive hypothesis. 
This concludes the case j = m+ 1. Now let j = 1, . . . , m, and suppose that 

ker<^ ("1 G i]+ i — ■ ■ ■ — kerip n G 4j+1 -i = kerip n G tj+1 = Kj+i 

In this case we have to show that for every I = ij-i + 1, . . . , ij, 

ker<^ n Gi = K 3 

and again we proceed by induction on I. Let us just consider the base case I = ij, the 
inductive step being similar to the one in the case j = m + 1. Let g G ker<p fl Gi j . 
Then g = gf ji for some e G Z and some h G G^+i. If e = then g G G^.+i 
and we conclude by inductive hypothesis that g G Kj+i- Now let us suppose 
e =/= 0. Arguing as before, the relation lattice contains an element of height ij and 
leading coefficient e. Thus e^' divides e. Let us denote by / the quotient. Then 

gGy+i = kf Gy+i, hence by inductive hypothesis gkj G ker^ n G^.+i = Kj+i, 
hence finally g G Kj. □ 

6. Some algorithms for lattices 

In this section we describe some algorithms that solve several problems related 
to lattices. We mainly work with matrices whose rows span lattices or subspaces 
in Q n . We say that a matrix is integral if it has integer entries. 

A basic algorithm that we use is the Smith normal form: given an m x n integral 
matrix A this algorithm finds an to x n integral matrix S, and integral unimodular 
square matrices P and Q with 

(1) S is in Smith normal form (this means that there is an r such that di — 
S(i, i) is positive for 1 < i < r, S has no other nonzero entries, and di 
divides di+\ for 1 < i < r), 

(2) S = PAQ. 

For details on this algorithm we refer to [TT], §8.3. One property that we note is 
the following. 

Lemma 6.1. Let qi,...,q n denote the rows of Q^ 1 . They form a basis ofTU 1 , and 
di is the smallest non-negative integer such that diqi lies in the span of the rows of 
A for 1 < i < r. 

We also need an algorithm that appears to be well known: in the computer 
algebra system Magma ([3]) it is implemented under the name of "saturation". 
However, we have not been able to find a reference for it in the literature. For this 
reason we sketch a solution here. Let A be an to x n-matrix with integer entries. 
Let V C Q ra be the Q-space spanned by the rows of A. The problem is to find a 
Z-basis for the lattice Z" l~l V. Without loss of generality we assume that the rows 
of A are linearly independent. The key observation is the following. Let B be an 
to x n integral matrix whose rows span V . Then its rows span Z ra n V if and only 
if the Smith normal form of B has diagonal entries that are all equal to 1. This 
follows from Lemma |6. II This yields the following algorithm. 

Algorithm 1 (Saturation). Input: an to x n integral matrix A with linearly 
independent rows. 

Output: an to x n integral matrix B whose rows span Z™ n V, where V G Q™ is 
the Q-space spanned by the rows of A. 

(1) Let S, P, Q be the output of the Smith normal form algorithm with input 
A. 

(2) Let 5" be the matrix obtained from S by setting the diagonal entries equal 
to 1. 
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(3) Return B = P~ 1 S'Q~ 1 . 

Algorithm 2 (Intersection of lattice and subspace). Input: innxn integral ma- 
trix A whose rows span the full dimensional lattice L in Q", and an m x n matrix 
B whose rows span an m-dimensional Q-subspace W of Q™. 

Output: annxn integral matrix whose rows span L, and whose first m rows span 
the lattice W n L. 

(1) Let ei, . . . , e„ and b\,...,b m denote the rows of A and B respectively. Write 
b% = Y^j=i Pij e ji an< i l et B' — (Pij)> i 1 necessary multiply the rows of B' 
by integers in order to get integral entries. 

(2) Let C be the output of Algorithm [1] with input B'. 

(3) Let S, P, Q be the output of the Smith normal form algorithm with input 
C. 

(4) return Q~ l A. 

Lemma 6.2. Algorithm^ is correct. 

Proof. The idea is to use the given basis of L as a basis of Q". Let ip : Q" — > Q™ be 
the corresponding isomorphism. So if v £ Q" then ip{v) is the vector that contains 
the coefficients of v with respect to the basis of L. So after the first step the rows 
of B' form a basis of tp(W). Of course tp(L) = Z" C Q™. 

So the rows of C form a basis of ip(W) O ip(L) = ip(W D L). Furthermore, the 
Smith normal form S of C has diagonal entries equal to 1. Therefore the rows of 
Q^ 1 form a basis of Z™ and the first m rows form a Z-basis of ip(WC\L). Note that 
for v £ Q™ we have = tp(v)A. Therefore the rows of Q _1 A form a basis of 

L, and the first m are a Z-basis oi W P\ L. □ 

Algorithm 3 (L-complements) . Input: an n x n integral matrix A whose rows 
span a full-dimensional lattice L in V = Q n ; and bases of subspaces V\ C y„_i C V. 
Output: an n x n integral matrix C with the following properties: 

- The rows of C span L. 

- The first s rows of C span Ln Vi (s = dim V\). 

- The first i rows of C span L n V n -\ {t = dim(V^_i)). 

(1) Execute Algorithm [2] with input A and a matrix whose rows span V n -i. 
Let Wi, . . . ,w n denote the rows of the output. 

(2) Let vt,...,v a be the given basis of V\, and write Vi = y] ■_ 1 ctijWj. Let 

(3) Execute Algorithm [2] with input the t x i-identity matrix, and A'. Let B 
denote the output. 

(4) Let C be the product of B and the t x n matrix whose rows are w\, . . . ,w t . 
Let C be the matrix obtained from C by appending w t +\, . . . , ui n . 

Lemma 6.3. Algorithm^ is correct. Let u\, . . . ,u n denote the rows of its out- 
put matrix B. Let Wi, W n -\ C Q™ be the subspaces spanned by m s +i, . . . , u n and 
Ut+i,...,u n respectively. Then W\ C W n —\ are a system of L-complements to 

Proof. After the first step w\, . . . , w n span L and wi, . . . ,w t span L n V n -i- Next 
we work in V n -\ using the basis Wi, . . . , Wt- We rewrite the basis elements of V\ 
with respect to this basis. We note that multiplying a v £ Q with the matrix with 
rows u>i, . . . , Wt is the reverse transformation. So the rows of C span L PI V^_i and 
the first s rows of C span Ln Vi- Therefore the output is correct. 

The last statement follows directly from the definition of L-complements. □ 
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Remark 6.4. The output of Algorithm [3] has one more useful property. Let 
W\,W n -i be as in Lemma IB31 The bases of the spaces W\ and W n -i that are 
produced by this algorithm are bases of L n W\ and L n W n -\ respectively. 

Algorithm 4 (integral relations). Input: An m x n-matrix A with rational coef- 
ficients. 

Output: an in x n integral matrix whose rows are a basis of the lattice 

m 

A = {(e ll ...,geZ m |^e i a i eZ"}, 

i=l 

where a\, . . . , a m are the rows of A. 

(1) Let M be the matrix obtained by appending the n x n-identity matrix at 
the bottom of A. 

(2) Let vi, . . . ,v m be a basis of the space {v € Q m+n \ vM = 0}. If necessary 
multiply each Vi by an integer to ensure that it has integral coefficients. 

(3) Let B be the output of the saturation algorithm (Algorithm [T]) applied to 
the matrix with the Vi as rows. 

(4) Output the rows of B with the last n coefficients deleted. 

Lemma 6.5. Algorithm^is correct. 

Proof. Note that the matrix M has rank n; therefore in Step 2. we find m linearly 
independent basis vectors. Now set 

A' = {e = (ei, . . . , e m+ „) 6 Z m+n \ eM = 0}. 

Then (ei, . . . , e m + n ) i— » (ei, . . . , e m ) is a bijection A — > A'. Now after Step 3. B is 
a basis of A'. We conclude that the output is a basis of A. □ 

7. The main algorithm 

Now we return to our initial problem. Let G C GL m (C) be a unipotent algebraic 
group defined over Q. Set V = Q m and let L be a full-dimensional lattice in V. 
The problem is to compute a finite set of generators of the group Gl = {g £ G(Q) 
g(L) = L}. 

We assume that G is defined as subset of GL m (C) by polynomial equations 
that have coefficients in Q. Then as a first step we find the Lie algebra of G 
as follows. First we compute a set S of generators for the radical of the ideal 
generated by the polynomials that define G (cf. pQ). Then we obtain the Lie 
algebra by differentiating the elements of S. 

The second step will be to compute a flag = Vq < V\ < ■ ■ ■ < V n = V of V for 
the action of 0. This is done by straightforward linear algebra: V\ is the space killed 
by all elements of 0; V2/V1 is the subspace of V/Vi that is killed by all elements of 
0, and so on. 

Now we have the input for our main algorithm which we now state. We use the 
notation of Sections [H |U 

Algorithm 5 (Main algorithm). Input: a non-zero finite dimensional vector space 
V over Q, a full-dimensional lattice L of V, a Lie subalgebra C fll(V) that is the 
Lie algebra of a unipotent algebraic group G, and a flag 

= V < Vi < • ■ • < V n = V 

of V with respect to the action of 0. 
Output: a T-sequence for Gl- 

(1) If n = 1 then return the empty set, else go to Step 2. 
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(2) Compute the derived vector space V*, the derived lattice L* , the derived 
flag = F * < VJ" < • • • < V*_ Y = V*. 

(3) Compute the kernel n and the image q of the derived action of jj on V* , 
together with the projection dir : $ — > q. 

(4) Apply the algorithm recursively to the vector space V*, the lattice L*, the 
Lie algebra q and the derived flag. Denote by gi, . . . , qk the result. 

(5) Compute a system W n -\ Q W\ of L-complements to V\ C V n -\, (Algo- 
rithm [3J the induced lattice Y and the induced error map e : End(l^) — > 
Lin(W„_i,Vi). 

(6) Compute a basis x%, . . . , xi of tlx, and set rij = exp(xi) for 1 < i < Z. 

(7) For 1 < i < k compute a preimage x of log(gj) through d7r and set g 9i = 
exp(x). 

(8) Compute the image W of n through e. 

(9) Compute a basis W = {wx, . . . , Wk} in Hermite normal form of the relation 
lattice of the elements e(g 9i ) +T + W in Un( ^ w 1 ,' Vl) for i = 1, . . . , k. 

(10) For each W{ in W do the following 

(a) Write Wi = (e^\ . . . , e^) and set 

3w>i — 3<?i ' ' ' ff<3/t • 

(b) Compute v Wi G W and 7 Wi £ T such that v Wi + r y Wi = e(g Wi ). 

(c) Compute the preimage n Wi of — v Wi through e : n — > Lin(W n _i, Vi). 

(d) Compute ^ = • exp(n TOi ). 

(11) Return gi, ... ,g fc ,rti, ... ,rt ; . 

We start by commenting on the computability of various steps. 

5. We note that the output of Algorithm [3] contains bases of L n W n -i and 
L n Vi. We uses these bases to represent an element of Lin(Wn-i,Vi) as 
an s x f-matrix (where s = dim(Vi), t = dim(Wn_i)). Then a Z-basis of T 
is the set of elementary s x i-matrices, which have one coefficient equal to 
1, and all other coefficients equal to 0. Computing e(a) for an a G End(V") 
is standard linear algebra. Indeed, for a G End(U) and w G W n -i write 
aw = V\ + Wi, where v\ G V± and w\ G W\. Then e(a)w = v±. 

6. Here we first compute a basis of the space e(n) C Lia(W n -i, Vi). Using 
Algorithm [2] we find a Z-basis of T n e(n). The inverse images of the basis 
elements under e are then a basis of (Proposition 14. 2|l . 

8. In Step 6. we already obtained a basis of W, and a basis 71,.. ■ ,7st of 
T such that 71, . . . , 7 m are a basis of W H T. Let A = si — m and let 
V> : Lin(W„_i, Vi) -> defined as follows. For 7 G Lin(W„_i, Vi) write 
7 = Z)i Cl >; then ^1(7) = (c m+ i, . . . , c st ). Then ^ is linear and 7 G T + W 
if and only if ^(7) G Z w . Set m — e(g qi ) for 1 < i < k. We want to 
compute a Z-basis of the lattice 

fc 

A = {(ei, . . . , e fc ) G Z fe I J2 e ^ ^ w ' + r l- 

i=l 

Now (ei, . . . , efc) lies in A if and only if J7 eiV( M i) € Z^. So we get a 
basis of A by applying Algorithm U] with input the matrix with rows ip(ui). 
Then we compute the Hermite normal form (cf. [TT], §8.1) of the basis of 
A obtained. 

The other steps are straightforward. We now prove the correctness of the algo- 
rithm. 
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Theorem 7.1 (Correctness). Let V be a non-zero finite dimensional vector space 
over Q, L a full- dimensional lattice of V , and G a unipotent algebraic subgroup of 
GL{V). Further, let q C qI(V) be the Lie algebra of G. Then Algorithm^ with 
input V , L and q and a flag 

= V < Vi < ■ ■ ■ < V n = V 

of V with respect to the action of g, returns a T -sequence for the subgroup Gl of 
GL(V). 

Proof. As a first thing we notice that, since g is the Lie algebra of G, as seen in 
Section 2] the given flag is also a flag of V with respect to the action G. Hence 
the hypothesis of Section [2] is satisfied, and we can consider all the constructions 
described there. As seen in Section n C g[(V) is the Lie algebra of N and 
q C flt(V*) is the Lie algebra of Q. Now we proceed by induction on the length n 
of the flag. 

If n = 1, then every vector of V is G-fixed, hence G is the trivial subgroup 
of GL(V) and Gl is the trivial subgroup of GL(V), hence the empty set is a T- 
sequence for Gl- 

Now suppose that n > 2. By the inductive hypothesis, qi, . . . , <j% is a T-sequence 
for the subgroup Ql* of GL(V*). By results in Section 01 the exponential map 
gives a bijection from n to N(Q) and, since by Proposition 14.11 the Lie algebra n is 
abelian, it is also a group morphism. Further, using Proposition 14.31 and Lemma 
12.51 it is easily seen that the image of through the exponential map is Nl . Hence 
ni, . . . , n s is a T-sequence for Nl , no matter in which order they are taken. As 
seen in Section[H we have exp(d7r(x)) = 7r(exp(a;)) for all x £ Q. So for i = 1, . . . , k 
we get 

(4) <g q ,) = q l - 

In particular, g qi G Gl* ■ By commutativity of the diagram in Corollary 1 2. A\ 

*(qi)=e(g qi )+r + W 

and, due to Proposition ^. 31 W is equal to the image W of N(Q) through e. Hence 
W is a basis in Hcrmitc normal form for the relation lattice of ^(qi), i — 1, . . . , k 
in Lin( ^ 1,Vl) . For 1 < i < k set 

hi = qi ■■■ q k k , 

where Wi = (e^\ . . . , ej^ ) is as in the algorithm. By Lemma 1 5. 11 the ordered set 
hi, . . . , hk is a T-sequence for the kernel of \&. Due to ((3]), we have that g Wi is an 
element of Gl* satisfying 

^(Sttu) = hi 

for 1 < i < k. Again by commutativity of the diagram in Corollarv l2.4[ we have that 
e(sftut) eT + W = T + W (hence Step HObl makes sense). Since exp(n Wi ) G N(Q), 
we have gi e Gl* and 

n{9i) = hi. 

Further, using Propositions 12.31 and 14.31 we obtain that 

e(gi) = e(g Wi ) + e(n Wi ) = v Wi - v Wi + j Wi G T 

hence by Lemma 12.51 we have that gi G Gl ■ Thus by Proposition 12.71 the ordered 
set QiNl, ■ ■ ■ , gk^L is a T-sequence for j^. Using Proposition 14.11 we get that 
N is central in G, hence Nl is central in Gl- Therefore we finally obtain that 
gi, . . . ,gk, rii, . . . , rii is a T-sequence for Gl- Q 
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Corollary 7.2. Let the notation be as in Theorem \7. 1\ Let TL denote the T- group 
generated by the output of Algorithm^ Then the Hirsh length ofTL is equal to the 
dimension ofG. Moreover, the Lie algebra of the radicable hull ofTL is isomorphic 
over Q to 0. 

Proof. We use the notation of the proof of Theorem l7.ll By induction, qi, . . . , q^ is a 
T-sequence for Ql* ■ This group has dimension equal to dim G — dim N = dim G — l. 
Therefore the T-sequence output by the algorithm is of length equal to dim G. 

Let hi,...,h r denote the T-sequence output by the algorithm. Due to jTtT| 
Chapter 6, the Lie algebra of the radicable hull of TC is isomorphic to the Lie 
algebra spanned by \og(hi). But the latter one is 0. □ 

Remark 7.3. The algorithm can be slightly modified in order to compute even a 
finite presentation of Gl- To show how it can be done, we need to introduce some 
notation first. 

Let Q be a T-group and g\, . . . , g n a T-sequence for Q. We say that a word w in 
gx, ■ ■ ■ ,g n is normal if it is of the form g^ 1 ■ ■ ■ g^ 1 for some £ Z. If this is the case, 
the depth of w is the minimum i such that is non-zero. It is an easy induction to 
prove that every element g of Q can be written in a unique way as a normal word w 
in g\, . . . , g n . Also, g € Gi = (gi, ■ . ■ , g n ) if and only if w has depth at least i. Since 
gx,...,g n is a T-sequence, for every 1 < i < j < n we have that [gi,gj] £ Qj+i> 
hence there exist a unique normal word W[ gug .] (of depth at least j + 1) such that 

[9i,9j] = 

It is well known that 

(gi, ■ ■ ■ ,9n | [gi, 9j] = W[ gii9j .] for I < i < j < n) 

is a finite presentation for G. We call it the standard presentation for Q with respect 
to gi, ...,g n . 

Recall that the algorithm as it stands computes a T-sequence g±, . . . , gk, ni, . . . , ni 
for Gl- Now we want to sketch how it can be modified in order to compute the 
standard presentation for Gl with respect to such a T-sequence. Of course it will 
be enough to show how to compute the normal words of the forms W[gi, gj ], u)[ gi ,nj] 
and W[ nitnj ] (for any suitable choice of the indexes i and j). The proof of Theorem 
17.11 shows that Nl is central in Gl ■ Hence it follows at once that 

W [9i,n 3 '] = W [n„n ] ] = 1. 

So the only hard part is to compute the words of the form u)\ gitg .y To this end 
we can suppose by inductive hypothesis that Step 4 of the algorithm, along with a 
T-sequence q\, . . . , qu for Ql*, produces also the standard presentation 

(qi, fe, qj] = w [q ^ qj] for 1 < i < j < k) 

of Ql* with respect to q±, . . . , q^- Recall that, using the notations of the proof of 
Theorem 1 7. 1[ hi, . . . , hk is a T-sequence for the kernel of Further, the proof even 
shows that we can effectively write its elements as normal words in qy , . . . , (of 
increasing depth). There are algorithms known for computing a standard presenta- 
tion of a subgroup of a T-group (cf. llj; an implementation of the algorithms for 
this purpose is available in the GAP4 package "polycyclic" , |7]). So we can compute 
the standard presentation 

(h 1 ,...,h k \ [hi, hj] = w [hi:hj ] for 1 < i < j < k) 

for the kernel of \& with respect to h\,...,hk- Since every Wr/ li ,/ l i is a normal 
word in hi,...,hk, we can evaluate it in gi,...,gk (that it to say, substituting 
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every hi with gi). In this way we obtain a normal word U[g i>g ■] in g\, . . . , g n . Since 
7r : Gl ™> K sends <?; in hi and has kernel Nl, we obtain that 

Since rii, . . . ,n s is a T-sequence for Nl, we can write n as a normal word v in 
ni,...,n s . This can be done effectively, being equivalent to write an element 
of a lattice in terms of a (ordered) basis. Since U[ gijg .jV is a normal word in 
gi, . . . , gu , ni , . . . , n;, it now follows easily that 

= 

8. Practical performance 

It is rather straightforward to see that the complexity of Algorithm is expo- 
nential in the length of the flag of V. Indeed, if the flag has maximal length, or in 
other words, dim V% = i, then dim V* = 2n — 2. So in the worst case the dimension 
of the ambient vector space is roughly doubled in each step of the recursion. 

It turns out that the dimension of the spaces Lin(W„_i, Vi) increases even faster 
(in the worst case by about a factor of 4 each step of the recursion). For this 
reason we avoid working with the entire space Lin(TVn-i, Vi). Instead we consider 
the associative algebra with one A C End(V) generated by the elements of the 
Lie algebra 0. Let U — e(A) C Lin(W„_i, Vi). Then e(G) C U, so we can work 
with the space U instead of Lin(W„_i, Vi). In fact we choose to work with a po- 
tentially somewhat bigger space, namely the subspace of Lin(W„_i, Vi) consisting 
of the matrices that have nonzero entries only in those positions for which there 
are elements in U that have nonzero entries in those positions. This space has the 
advantage that its intersection with the lattice T is easily computed. Furthermore, 
in practice it is only slightly bigger than U. 

We have implemented the algorithm in the language of GAP4 ([8 ). We use two 
series of Lie algebras in gI n (Q) to generate test inputs to the algorithm. The terms 
of the first series are denoted 0„ , which consists of the elements 

%i = for 1 < i < n — 1 

n-1 

(here e^j is the n x rt-matrix with a 1 on position and zeros elsewhere). The 
only nonzero commutators are [xj,a; n ] = Xi + \ for 1 < i < n — 2. So 0„ is of 
dimension n and of nilpotency class n — 1. 

The second series of Lie algebras is denoted f)„, which consists of the elements 

n-1 

Vi = ^ «e M+ i 
i=l 
n—k 

Vk = e i,i+k for 2 < fe < n — 1. 

i=l 

Here the only nonzero commutators are [j/i, yk] = —ky^+i for 2 < k < n — 2. So t) n 
is of dimension n — 1 and of nipotency class n — 2. In fact, as abstract Lie algebras, 
f)„ = n -i- We note that both Lie algebras have a flag of maximal length. 

In Table[T]we list the runing times((j) of the algorithm with input the Lie algebras 
9 n and h n , for n — 6,7, 8, 9. In all cases for the lattice L we have taken Z" C Q™. 

From Table [T] we see that the algorithm is efficient enough to tackle nontrivial 
examples. However, the running times do confirm the analysis above that the 



The computations were done on a 2GHz processor with 1GB of memory for GAP. 
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n 


time 9 n 


time f)„ 


6 


0.7 


0.4 


7 


3 


3 


8 


24 


16 


9 


204 


133 



Table 1 . Time (in seconds) for the main algorithm with input 3 
and t) n . 



complextity of the algorithm is exponential. Also we see that the running time is 
essentially determined by the length of the flag, as fj n = S n -i, but the algorithm 
needs markedly longer for f)„ than for 0„_i. 
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